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A system with a topological insulator slab sandwiched by two one-dimcnsional photonic crystals as 
a defect is investigated. A giant topological magnetoelectric effect is proposed which is characterized 
by two resonant peaks with half-unit transmissions and rich polarization behaviors in the photonic 
band gap. We also predict a giant optical Hall effect in our system, in which the transverse shift of 
the transmitted light could be thousand times larger than the wavelength. These effects are purely 
from the 0-domain walls at the topological insulator interfaces and could be observed easily. These 
effects can also be used as the sensitive detection of basic physical constants. 
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Introduction. The modification to the Maxwell La- 
grangian of classical electromagnetism by a term, Ca — 
E • B, is known as "axion electrodynamics" in field 
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theory literature where a = /he is the fine structure 
constant, and 8 is known as the axion field. Novel effects 
arise when the axion field is a function of space-time 
0]. Recently it is argued that Ca can be realized as the 
low-energy physics of three-dimensional (3D) topological 
insulators (TI) which have been theoretically pre- 
dicted and/or experimentally observed in various binary 
[i-d] and ternary [i-|T2| compounds. In this context, the 
Lagrangian Ca describes the topological part for the elec- 
tromagnetic response of 3D TI with modified constitution 
equations D = eE - a(e/7r)B, H = B/^ + Q;(e/7r)E, 
where e and /i are respectively the material permittivity 
and permeability. In the above, = or (±)7r is viewed 
as a phenomenological parameter which categorizes all 
time-reversal invariant insulators into Z2 trivial or non- 
trivial class. Since the term can be written as a total 
derivative in the Lagrangian, it causes nontrivial effects 
only at the interfaces, which separate TI and normal in- 
sulators and are also called as 8-domain walls (DWs). 
These effects on 8-DWs are known as the topological 
magnetoelectric effects (TME), in which an electric field 
can induce a magnetic polarization and a magnetic field 
can induce an electric polarization. To eliminate the sign 
ambiguity of for Z2 -nontrivial insulators, a T-breaking 
surface magnetic field is applied to determine the wind- 
ing direction of through the DW, so that 8 = tt when 
the magnetization is in the interface normal direction and 
— TT when the magnetization is inverted. 



As a direct consequence of Z2 nontrivialit 
tentions have been focused on the TME 13 



V, g reat at- 
16|, among 

which one jl3l | abandons the toy semi-infinite geometry 



* liuqin@mail.sim.ac.cn 
txyjiang@mail.sim.ac.cn 



9r 



I J - 



□ 



glass 



TI 



FIG. 1: Schematic show of our model with a TI slab as defect 
sandwiched by two ID PhCs. 



and seriously considers the realistic thin-film structure of 
TI [l3| as well as the weak exchange- coupling to the sur- 
face magnetization. In this connection, optical methods 
take special advantages in measuring the TME experi- 
mentally [3 [3 J e-ff- by measuring the Kerr or Faraday 
rotation angles, 6'(B), in the zero B-field limit. However, 
TME is very weak generally since 6'(i? — >■ 0+) cx; a and 
has not been observed experimentally so far, hence, find- 
ing larger physical effects to manifest Z2 nontriviality is 
very essential. 

Motivated by these facts, in this letter, we study the 
transmission characteristics of a system with a 3D TI slab 
sandwiched by two ID photonic crystals (PhCs) as a de- 
fect using the standard transfer-matrix method (TMM) 
[iol ]. It is found that, there are two resonant peaks in 
the photonic band gap (PBG) with half-unit transmis- 
sion each, which is in sharp contrast with a single uni- 
modular resonant peak of a normal defect. For normal 
incidence with linear polarization, the transmitted waves 
at the two resonant peaks are left and right circular po- 
larized respectively, and at the central frequency between 
the two peaks, the polarization of the transmitted wave 
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is linear but with the Faraday rotation angle about 7r/2. 
These diverse phenomena play the role of "giant TME" 
and can be optically detected at Terahertz (THz) for real 
TI materials. Moreover, with small incident angle, a gi- 
ant transverse shift, which could be thousand times larger 
than wavelength, is predicted for the resonant transmit- 
ted light, which is termed as "giant optical Hall effect" 
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Our system, as shown in Fig[Tl is made of a FM-TI- 
FM heterostructure and two symmetric N-cell ID PhCs 
at both sides. Each cell consists of a glass and a vacuum 
layer. The surface FM magnetization in 2-direction gen- 
erally could be parallel and antiparallel in which the 
0-DWs at the left and right surfaces of the TI are re- 
spectively (0,7r)-(-7r,0) and (0, 7r)-(7r, 0). The FM layer 
opens an energy gap of the surface states of TI, which is 
typically Eg = 10 meV [^j. Since the effective theory Ca 
applies only in the low frequency limit, ui <^ Eg/h, which 
implies that the working frequency of interest is within 
THz range, to detect the topological phenomena, the re- 
fractive indexes and the thickness for the binary PhC 
cell and the TI slab are taken respectively as Ug = 1.5, 
riv = 1.0, riTi = 10, and d = dg = dy = dn =0.1 mm. 
Under this choice, the central frequency of the first PEG 
is around 0.6 THz, and the finite size effect of the TI 
defect can be neglected [23 1 . 

Giant topological magnetoelectric effect. As the start- 
ing point, we first consider the parallel magnetization 
case where = — ©fl = tt. The transmission spec- 
trum with PhC cell number = 10 for a normally in- 
cident linear plane wave with Eq along x-axis is shown 
in Figl2][a). Surprisingly, two topological-defect resonant 
peaks are observed with half-unit transmission each. We 
have also checked other cell numbers and the evolution 
of peaks with N is shown in Figl^l^b). It is found that 
when N < 5, there is a single almost- unimodular reso- 
nant peak in PEG, with increasing N, the single peak 
splits into two, which turn to be so sharp that they are 
well-separated when N > 7. Furthermore, at N = 10, 
the polarization of the transmitted light changes from 
right (wi) to left elliptic (^5) with the frequency increas- 
ing as in FiglSJc). There arc three frequencies at which 
the polarizations are of particular interests. At two reso- 
nant frequencies W2 and UJ4, the transmitted lights are of 
right and left circular polarization respectively, while at 
the central frequency W3 between them, the transmitted 
light is linearly polarized again but with a giant Faraday 
rotation angle about 7r/2. 

To understand these anomalous phenomena of topolog- 
ical defect in PhC, we set up a simple Fabry-Perot (FP) 
cavity model for our system. As shown in Figll] if we take 
the areas in the two dashed boxes (which include the 0- 
DWs at the surfaces of the TI slab and PhC) separately 
as left and right effective mirrors, then the system can 
be treated as a typical FP cavity. Before taking deeper 
insight into the TI FP cavity, we first review the theory 
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FIG. 2: (a) Transmission h/Ii and reflection Ir/h spectra 
with incidence as a linearly-polarized plane wave, (b) Evolu- 
tion of transmission peaks with cell number A''. The solid lines 
are numerical results, while the dash-dots are from our the- 
oretical resonant conditions, (c) Angle (j) (blue line) between 
the long-axis of polarization and a;-axis, and the intensity ra- 
tio (green line) between the short- and long-axis I = h/Ia 
as a function of frequency. The insets are schematic polar- 
ization of transmitted wave at five typical frequencies, where 
the long-axis is shown in red. 



of a conventional one [20| . For a FP cavity with conven- 
tional material, the transmission coefficient for the FP 



cavity is well-known as tpp = iit2e - /(I — rir2e ' ), 
where ti, ri and t2, ^2 are respectively the left-to-right 
transmission and reflection coefficients of the two mir- 
rors, and kzd is the one-way optical path through the 
cavity. This FP transmission coefficient is usually pic- 
tured as the sum of all orders back-forward scattering 
between the two mirrors, tpp = tit2e^^''^\l-\-rir2e'^''^''^ + 
{rir2e'^^^''^)'^ + ■••], from which we see that the resonant 
condition + 4'r2 + '2kz(Ld)d = 2m7r. m e N, is satisfied 
simultaneously for both TM = (1 0)'^ and TE = (0 1)"^ 
modes, where ri_2 = |''i,2| ^^p(J'4>ri,r2)- Therefore there is 
one single resonant peak in transmission spectrum with 
unit transmission Tpp = |tFPp- 

However for a TI FP cavity, the above theory needs 
to be modified due to the existence of 6-DWs. In previ- 
ous studies [i, 13, 18 1, nonzero Kerr or Faraday rotations 
at 8-DWs are predicted, where waves of different linear 
polarization couple to each other. This fact tells us that 
the TM and TE modes are no longer good bases in the 
presence of 0-DWs, and in these bases, the reflection 
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(transmission) properties are described by a 2 by 2 ma- 

trix, termed as "reflection (transmission) matrix" r [t). 
Therefore, to obtain the resonant condition following the 
same discussions as the above, the key question is: what 
are the good bases for systems with 0-DWs? This prob- 
lem is resolved by finding the form of 8-DW reflection 
(transmission) matrix as: 



r (0 = 



_ f ^r(t) B,.(t) 



(1) 



where A^i^t) and Bj.{t) are real, and the off-diagonal term 
generated by the 9-DW, i? oc a, is a small value. 
The good bases of such matrix are nothing else but 
a± ~ {1 ^ i)"^/\/2, which are respectively the right and 
left circular polarized waves. Physically this means, if 
we normally incident one eigenmode of r , say cr+ , then 
the reflected wave is still a+ but multiplied by a reflect- 
ing coefficient row = \'r'Dw\&'>^pii4>Dw), where Iruvyl = 
(A^ + -8^)1/2 ^ \Ar\ and (j)DW = arctan(-5r/Ar) is a 
topological phase generated by 6-DW. For our TI FP 
cavity where each effective mirror includes not only a 
0-DW but also a PhC, the essential tricky is that the re- 
flection (transmission) matrix of such mirrors still has 
the form of Eq.(IT]), but with Aj. = \Ar\exp{i(pr) be- 
ing a complex number and 4>dw = arctan[5R(— B^/^r)]- 
Therefore, a± are still good bases and any normal in- 
cident wave could be expanded as Eq = {E+a+ + 
E-a-)/V^, where E± = ± iEy, so the discussions 
for a conventional FP cavity are valid separately for 
each basis of (7±, and the transmission of our system 
■ rp — 1 V (1--R)'' 

(l--R)^-|-4_Rsin'' I k^d+ct>r+s i-j ^ 



where i? = + l-Brp. Obviously, the theoretical res- 
onant conditions for a± modes are: 

2k^{u)d + 2(j)r ± {(t)DWi + ^uws) = 2m7r, m G N. (2) 

The phases ±{(I)dwi + (f>DW2) generated by two 9-DWs 
lead to two resonant frequencies 0^2,4 for a± modes. In 
FiglUb), we plot the resonant frequencies obtained by 
Eq.® for different N in dash-dots, which is compared 
with those obtained by numerical methods shown in 
solids, and we sec that they agree with each other very 
well. Two essential properties are concluded. First, the 
PhCs are crucial to observe the double peaks due to their 
narrowing down peakwidth dramatically. Second, the 
frequency difference between the two peaks, Suj = UJ4—UJ2, 
is almost a constant as the change of N, although the 
two resonant peaks are unsolvable in spectra for small 
A^. Actually, even without PhCs, the resonant condition 
Eq.® still gives two resonant frequencies for a bare TI 
slab, so the two-peaks spectrum is an intrinsic property 
generated by 0-DWs. For our system, it is found that 
i5a;/w ~ , a"^^2-i Ct- Since experimentally we can mea- 
sure the frequency almost exactly by interferometers, Slj 



is a new way not only to manifest topological nontriv- 
iality but also to measure the basic physical constant 
a = /he. 

In Figl^Kc), a rich behavior of polarization is shown 
versus frequency, which can also be understood in the TI 
FP cavity model. At resonant frequency, W2 or a;4, only 
one eigenmode is totally transmitted while the other is 
totally reflected, so the total transmitted light is circular 
polarized. While at wa the transmission amplitudes of 
both a± modes are small but the same, so that the total 
transmitted wave is linearly polarized again, but with the 
Faraday rotation angle about 7r/2 (or 37r/2). 

In contrast, if the TI slab has antiparallel magnetiza- 
tion at two interfaces with Ql = 0_r = tt, the TI FP 
cavity goes back to a conventional one with a single uni- 
modular resonant peak. This is because in this case, 
A-r 



Ar2 while Bj-i 



'Br2, SO that 



'JDW2 



and Eq.® recovers that of the conventional FP cavity. 
Unlike the parallel magnetization case where the Faraday 
rotation is doubled when light passes two 8-DWs from 
left to right, for antiparallel magnetization, the opposite 
half-quantized Hall conductances carried by the two sur- 
faces cancel each other, hence there is no net Faraday 
rotation when out of the TI defect. Vanishing of TME 
in antiparallel case is another demonstration of pure Z2- 
nontriviality effect of our system. For experimentalists, 
antiparallel case is also a way to eliminate systematic 
errors. 

Giant optical Hall effect (OHE) of TI defect. In elec- 
tronic systems, the anomalous velocity (k x $7^) origi- 
nating from Berry curvature ilt- l 24|.l25i in various Hall 
effects has been well studied 12614281 . In parallel On- 
oda et al. construct in Rcfs. [2l|, |22| a theory for the 
propagation of an optical wave packet in slowly- varying 
photonic systems, such as chirped PhC. A set of general 
equations of motion (EOM) is derived in the adiabatic 
limit, from which the anomalous velocity is proved to be 
consistent with a more basic physical requirement, i.e., 
the total angular momentum conservation (TAMC). The 
nonvanishing anomalous velocity predicts many interest- 
ing dynamical phenomena, one of them is the OHE which 
could be measured by a transverse shift of the transmit- 
ted light beam at the interface between different media as 
shown in Fig|3Ka). Strictly speaking, the EOM is not ap- 
plicable for sharp change at interface due to the failure of 
the adiabatic condition, however the correct beam-shift 
value can be obtained by the TAMC law. Generally, the 
transverse shift in OHE is only a fraction of the wave- 
length. A kind of giant OHE, where the transverse shift 
is thousand times of the wavelength (or 10^ atomic lattice 
constant), is predicted at the interface between a PhC 
and vacuum when a small PEG A is opened at the orig- 
inal Dirac point by spatial-inversion symmetry breaking 
[21I . The mechanism of the giant OHE is understood 
by TAMC such that a very strong pseudo-spin of Bloch 
function in PhC needs to be compensated by the orbital 
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angular momentum in vacuum in the form of large trans- 
verse shift. 

In our system, we report a completely different mech- 
anism to generate the giant OHE. Before getting into 
the details, let's first go over the calculations of the 
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OHE transverse shift at interface by TAMC 
For a system with rotational symmetry around the z- 
axis, the ^-component of the total angular momentum 
with both orbital and "spin" (including pseudo-spin of 
the Bloch function) degree of freedoms, = [fc x kc -I- 
{zc\<Jy\zi}){\<iclkc)\z, is conserved between incident and 
transmitted lights, j\ = jj. Here and kc are re- 
spectively the central position and wavevector, jzc) is a 
two-component spinor and is the Pauli matrix. The 
TAMC yields the transverse shift for the central position 
of transmitted light beam as 21 , 22] : 



1 

kl sin ( 



) cos 



(4k3Ni)cOS0/] ,(3) 



where ^/(t) is the angle between the z-axis and kc 
of the incident (transmitted) light. From Eq.(|21), the gi- 
ant OHE at PhC- vacuum interface 21], |22| can be ex- 
plained by the vanishing of the large PhC pseudo-spin, 
(Zclcal-^c) ^ l/^^i which diverges at the limit A — > 0. 
For our system, if we choose a linearly-polarized nearly- 
normal incident beam with frequency at the resonant 
one, W2 or 0^4, the transmitted light (with half-unit trans- 
mission) is right or left circularly polarized, which means 
that |[(2j|f73|zJ)cos6iT- (2c^|a3|4)cos6'/]| ~ |±1-0| = 
1 is almost a constant at small incident angles. How- 
ever, (5?/ J diverges since the denominator is approxi- 
mately l/{k^9i) — 00 in the limit 0/ — > 0. 

Using the TMM, the transmission and reflection fields 
under tilt incidence are also obtained, and Sy]^ as a func- 
tion of 9i, calculated from Eq.Q, is plotted in FiglJKb) 
where its divergence at small incident angles is clearly 
seen. 

Finally, it is interesting to discuss the physical limit 
for observing giant TME and OHE in our system. From 
the peakwidth of the resonant modes in Figl2ja), we note 
that the modern high quality lasers and interferometers, 
whose frequency deviation is in order of 10~^, can eas- 
ily catch the exact resonant frequencies 0^2(4) of our sys- 
tem, therefore giant TME can be observed. For giant 
OHE, the limit is not from the frequency regime, but 
mainly from the k regime. To observe the beam shift, 
the incident beam must have a finite- width Di,, which 
is composed of different wavevectors. Suppose the inci- 
dent beam is a Gaussian one, the distribution width of 
the wavevectors is Ak ~ 1/-Df,, which determines the 
minimum-achievable incident angle as 9J""' ~ Afc/fcc. 
Therefore the transverse shift has the magnitude of the 
beam width 5y ^ l/{k^9j) ^ Dt,, which could be nearly 
hundred or thousand times larger than the wavelength in 
the laboratory conditions for THz waves. Another way 
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FIG. 3: (a) Schematic show of OHE with transverse shift of 
transmitted beam, (b) Transverse shift versus incident angle 
at frequency (jJ2 for our system, (c) The spatial distribution 
intensity I{x, y) and the energy current direction of transmit- 
ted light when the incidence is a linearly-polarized Gaussian 
beam with width Di, — lOOA and frequency as ljJ2 in FigO 



to directly observe giant OHE is to detect the transmit- 
ted beam of a normally incident Gaussian beam to our 
system. Since the nearly-exact-normal fc-components of 
Gaussian beam have large OHE shift, we expect that the 
transmitted "beam" is a "light ring" whose center is dark 
and the maximum should shift to the position nearly Dt 
from the beam center as discussed above. Furthermore, 
since OHE is a transverse shift relative to incident in- 
terface, the energy current of transmitted beam should 
have the vortex form. Based on strict Green's function 
method, we have calculated the transmitted beam of an 
incident Gaussian beam with width as = lOOA = 5 
cm and the results are shown in Fig[3](c), which agrees 
with our expectation very well. Our discussion of both 
giant TME and OHE can also be applied to the reflected 
waves. At last, we note that the giant TME and OHE can 
be utilized not only as manifestation of topological non- 
triviality, but also as sensitive detecting methods with 
promising potentials since the simplicity of our system. 
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